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Abstract 

We study the diffusive motion of low-energy normal quasiparticles along the core of a 
single vortex in a dirty, type-II, s-wave superconductor. The physics of this system is argued 
to be described by a one-dimensional supersymmetric nonlinear a model, which differs from 
the (J models known for disordered metallic wires. For an isolated vortex and quasiparticle 
energies less than the Thouless energy -Bxh, we recover the spectral correlations that are 
predicted by random matrix theory for the universality class C. We then consider the 
transport problem of transmission of quasiparticles through a vortex connected to particle 
reservoirs at both ends. The transmittance at zero energy exhibits a weak localization 
correction reminiscent of quasi-one-dimensional metallic systems with symmetry index /3 = 
1. Weak localization disappears with increasing energy over a scale set by Bth- This 
crossover should be observable in measurements of the longitudinal heat conductivity of an 
ensemble of vortices under mesoscopic conditions. In the regime of strong localization, the 
localization length is shown to decrease by a factor of 8 as the quasiparticle energy goes to 
zero. 

PACS numbers: 74.40.+k, 74.25.Fy, 72.80.Ng, 73.23.-b 

1 Introduction 

Disordered single-electron systems fall into several universality classes distinguished by their 
symmetries under time reversal and rotation of the electron spin. In mesoscopic physics three 
such classes are well known, namely those anticipated by Dyson's 1962 classification ||] of many- 
body systems with complex interactions. They are called the classes with orthogonal, unitary, 
or symplectic symmetry. In the ergodic regime, where the relevant time scale is much longer 
than the time for an electron to diffuse across the system, each universality class is described 
by a Wigner-Dyson random matrix ensemble with symmetry index /3 = 1, 2, 4, respectively. 

In the present article we consider a fascinating variant of the theme of disordered electron 
systems: a single vortex in a disordered type-II conventional s-wave superconductor. A simple 
picture would model the vortex as a tube-shaped metallic region inside the superconductor. 
Indeed, a vortex supports quasiparticle excitations, which have normal characteristics and are 
bound to the vortex core by the pairing gap of the superconductor. It is then natural to ask 
how the quasiparticles of a disordered vortex fit into the traditional classification scheme. 
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The first and important message is that the low-energy quasiparticles exceed the traditional 
scheme and form a separate universality class. This assertion can be motivated as follows. 
Dyson's classification assumes the absence of symmetries of the particle-hole type. When this 
assumption is relaxed, seven more universality classes enter. Three of these are realized for a 
massless Dirac particle moving in a random gauge field Q , and another four have recently been 
identified Q in metallic systems in contact with a superconductor. Common to all of these 
classes is an invariance of the energy (or Dirac) spectrum under inversion E i— > —E. In the limit 
£ — > 00, the inversion symmetry becomes irrelevant and conventional Wigner-Dyson statistics 
is recovered, but for small enough E novel features appear. In particular, the disorder average 
of the density of states is nonstationary and yet universal. 

A vortex is an example of a metallic system with particle-hole symmetry. (Note that this 
particle-hole symmetry is not a dynamical symmetry but arises as a result of doubling the single- 
particle Hilbert space a la Bogoliubov-deGennes, in order to incorporate pairing correlations 
within the formalism of first quantization.) The magnetic field admitted by the vortex breaks 
time reversal invariance whereas, if magnetic impurities are absent, the spin of a quasiparticle 
is conserved. According to the classification scheme of [|| , this means that a disordered vortex 
belongs to symmetry class C. Moreover, the supercurrent circulating around a vortex causes 
the phase shift due to Andreev scattering to vanish on average, so that by the reasoning of |^ 
the Gaussian random matrix ensemble of type C is expected to be appropriate at low excitation 
energies. 

These predictions have recently been applied to dissipative vortex motion in a quick pro- 
posal (for pancake vortices) by Skvortsov and Feigel'man Q. The basic picture is as follows. 
Dissipation is due to the transfer of energy from the collective motion to the internal degrees 
of freedom. In the case of a vortex the relevant internal degrees of freedom are those of the 
quasiparticles. The rate of energy transfer is largely determined by the statistical properties 
of the quasiparticle wavefunctions and energies. Therefore, a modification of their statistics is 
expected to modify the rate of dissipation. Indeed, Skvortsov and Feigel'man argued that the 
level statistics of type C leads to a dissipative flux-flow conductivity which is weakly anomalous. 

Although random matrix theory in its universal form is a powerful and successful concept, it 
is limited to the ergodic regime of times larger than the diffusion time or, equivalently, energies 
smaller than the Thouless energy. At short times, or high energies, the diffusive dynamics of 
a disordered electron system causes deviations from random matrix statistics. This is true for 
ordinary metals and must equally be true for the quasiparticles of a vortex. Clearly, a complete 
theory of mesoscopic phenomena needs to encompass both regimes and make the crossover 
between them theoretically accessible. For the case of ordinary metals, such a theory exists 
and is in good shape. In its most elegant formulation, it is given by the diffusive nonlinear a 
model pioneered by Wegner |^ and perfected by Efetov Q. This is a nonlinear field theory of 
interacting diffusion modes, which has random matrix theory for its universal low-energy limit. 

Is it possible to establish a similar description for the quasiparticles of a disordered vortex? 
As follows from the above discussion, such a description is of relevance to the microscopic theory 
of dissipative vortex dynamics. It should also be directly testable by experiment, by measuring 
the longitudinal heat conductivity to probe the transport properties of the quasiparticles at low 
energies or temperatures. (Note that electrical conductance measurements are ruled out because 
the vortex is short cut by the superconducting condensate.) 

It is clear that the construction of a field theory for the disordered vortex, encompassing 
both the diffusive quasiparticle dynamics and its ergodic (or random matrix) long time limit, 
is a more intricate problem than deriving the nonlinear a model for disordered metals. This is 
because diffusion in a disordered superconductor is already a nontrivial nonlinear process. The 
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relevant kinetic equation was derived by Usadel , building on earlier work by Eilenberger g| . 
For the special geometry of a single vortex, the solution of the Usadel equation was obtained in 
Refs. ijl^ to predict the crude structure of the local density of states of the vortex and its 
thermodynamic properties. 

The Usadel equation is well established as a tool in the research area of disordered super- 
conductors and the proximity effect. However, because of its quasiclassical nature, it leaves out 
mesoscopic effects that are due to quantum coherence, such as weak localization and mesoscopic 
fluctuation phenomena. The question then is whether one can improve on the quasiclassical 
Usadel equation and incorporate the relevant quantum effects which, in the long time limit, give 
rise to random matrix features. The answer is yes. In an important recent letter, Altland et al. 
[ p] pointed out that the Usadel equation can be interpreted as the saddle point equation for a 
certain supersymmetric field theory of the type of a nonlinear a model. (The general connection 
between quasiclassical kinetics and nonlinear a models was first understood by Muzykantskii 
and Khmel'nitskii The formalism of Altland et al. also shows how to go further and 

include quantum interference effects, by computing the fluctuations around the saddle point 
solution. In the present paper this formalism will be applied to the case of a disordered vortex. 

A brief account of the contents of this paper is as follows. Section ^ defines the model and 
sets up the basic field theoretic formalism, the main step being the introduction of a composite 
field Q by means of a Hubbard-Stratonovich transformation. In Section ^ two saddle point 
approximations are made in sequence. The first of these restricts Q to a nonlinear field space 
obeying the local constraint Q^{x) = 1. As usual, this approximation is justified for diffusive 
systems, where the disorder is at the same time (i) sufficiently strong to cause rapid momentum 
relaxation (thereby leading to diffusive motion in coordinate space) and (ii) weak enough in 
order for the system to retain metallic characteristics. 

The second saddle point approximation is made in the course of solving for the dynamics of 
the Q field inside the nonlinear space = I. To that end, we expand the full Lagrangian of 
the theory to obtain the relevant "low-energy" Lagrangian, containing in addition to the usual 
diffusive (or kinetic) term a linear coupling to the superconducting order parameter, and the 
minimal coupling to the magnetic field present in the vortex. Such an expansion is valid in the 
so-called "dirty limit", where the coherence length (*-e- the length scale of variation of the 
pairing field) exceeds the elastic mean-free path The equation of motion obtained by varying 
the resulting action functional turns out to be the Usadel equation. Its solution depends only on 
the scaling variable D^^'^r, where r is the distance from the axis of the vortex, Ao is the 
magnitude of the pairing gap in the bulk, and D is the diffusion constant of the normal metal. 
This result follows on dimensional grounds and is in agreement with the well-known fact ||l3|| 
that the coherence length gets reduced by a factor {i/£,oY^'^ in the dirty limit. 

Section || explores the consequences of our key observation, which is that the solution of 
the Usadel equation is degenerate. More precisely, the field theory for Q has an exact global 
symmetry broken only by the energy distance from the Fermi level. (This symmetry is not 
present in traditional work based on the Usadel equation. We believe that its existence can 
readily be seen only when the field theoretic formalism is used.) Letting the symmetry group act 
on the solution of the Usadel equation we get a degenerate manifold of solutions. The "radius" 
of the manifold is maximal on the axis of the vortex and shrinks to zero with increasing distance, 
over a scale given by [D/ /S.^Y^'^ ~ [S.o^-Y^'^- Thus, the symmetry becomes ineffective as we move 
away from the vortex to the asymptotic region. In this way the symmetry, while being global or 
constant in space, manages to perform the curious trick of having an effect that is localized in 
space! 

For an isolated vortex of length L, the statistics of the quasiparticle levels at energies less 
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than the Thouless energy i?Th ~ D /L^ is determined by the manifold swept out by the symmetry 
group when acting on the solution of Usadel's equation. By integrating over it, we recover for 
the vortex the random matrix prediction for the level statistics of systems in universality class C. 
For higher energies {E > -Etii), or fixed energy and increasing system size L, long wave length 
fluctuations inside the manifold become important. These are governed by a one-dimensional 
field theory of the nonlinear a model type. We identify the field space as a Riemannian symmetric 
superspace of type DIII|CI. The field is called q. 

In Section ^ we turn to the investigation of transport properties of the quasiparticle exci- 
tations of the disordered vortex. We imagine the two ends of the vortex (or of an ensemble of 
vortices, to facilitate experiment) to be connected to infinite particle reservoirs and study the 
problem of transmission of quasiparticles along the axis of the vortex. We work in a mesoscopic 
regime, requiring low energy and temperature. This problem is the natural analog of calculating 
the disorder average of the conductance of a mesoscopic metallic wire. The specific observable 
we study is a correlation function for the conserved probability current carried by quasiparticles 
of energy E. We call this observable the "transmittance" , t{E). The coupling to the reservoirs 
is modeled by open boundary conditions on the field q in the standard way. 

We then carry out the usual field theoretic perturbation expansion up to one-loop order. To 
deal with the symmetry breaking due to finite energy E, which causes crossover between uni- 
versality classes, we introduce a parametrization of the fields akin to that used by Altland, lida 
and Efetov in their supersymmetric approach to the crossover problem for small disordered 
electron systems with orthogonal and unitary symmetry. The result we obtain is 



where ^loc has the meaning of a localization length, and e is the energy measured in units of 
(roughly speaking) the Thouless energy. For quasiparticles of zero energy, the result simplifies to 
r = {^loc/L) — (1/3) + which is similar to the conductance of quasi-one-dimensional metallic 
systems in the Wigner-Dyson universality class of orthogonal symmetry (/3 = 1). The second 
term in the expression for r is called the weak localization correction. It disappears in the limit 
e 3> 1, where our field theory reduces to Efetov's nonlinear cr model with symmetry index f3 — 2. 

Finally, in Section ^ we tackle the problem of exponential (or strong) localization, which 
emerges at lengths L/^ioc > 1- The natural tool for calculating the transmittance in this non- 
perturbative regime is the "quantum Hamiltonian" of the one-dimensional functional integral. 
By a standard transfer matrix argument, the exponential decay of the transmittance is deter- 
mined by the lowest nonzero eigenvalue of the quantum Hamiltonian. Calculating it for an 
arbitrary energy e is a problem beyond the scope of this paper. However, in the two limits e = 
and e S> 1 the quantum Hamiltonian is simply a Laplacian on a symmetric superspace and has 
a large enough group symmetry to be tractable. By computing the low lying spectrum of the 
Laplacian we find that the localization length for e = is eight times smaller than for e 3> 1. 

2 Model and disorder averaged Green functions 

Our starting point are the Bogoliubov-deGennes (BdG) equations for the stationary quasiparticle 
excitations [p^ of a superconductor in (BCS) mean field approximation: 




Hqu + Av = Eu , 
A*u - H^v = Ev . 
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The single-particle Hamiltonian Hq is expressed by 



Ho = i-iV - eAf/2m + V{x) - ^jl , 

where /j, is the chemical potential (Fermi level) and A is the magnetic vector potential. The 
functions u{x) and v{x) describe the particle and hole components of a quasiparticle, and are 
coupled by the pairing potential A(x). The potential energy V{x) is taken to be a Gaussian 
white noise random potential with zero mean and variance 

{V{x)V{y)) = vH{x - y) = Six - y) , 

where ly is the density of states and r is the elastic mean free time. It is convenient to assemble 
Hq and A into a single operator H, here called the BdG "Hamiltonian" for short: 



H^, Ho 



A* -H^ 



which acts on the tensor product of with the Hilbert space of square-integrable functions 
on the coordinate space, S, of the superconductor. The factor will be referred to as the 
particle-hole (ph) space. 

The BdG Hamiltonian has the following symmetry property: 



H = -CH^C^^ with C 



1 

-1 



where T means transposition in both particle-hole and coordinate space. As a result of this 
symmetry, the eigenvalues of H occur in pairs with opposite sign, and we have 



G{z) := {H - z)-^ = {-CH^C-^ - z)-^ 

= -C{{H + z)-^fc-^^~CG^{-z)C-\ 

so that the advanced and retarded Green functions are related by 

Gaa,{x,x';E) = - ^CatiCa'f}'G^,fj{x' ,X;-E) . 

To compute disorder averaged Green functions, we will use Efetov's supersymmetry method 
in its recent adaptation by Altland et al. |pl| to include the pairing correlations of a 
superconductor. In this method, disorder averages of Green functions are obtained from a 
generating functional 

Z[j] = ( y 'Di)Vij exp j (fx {ii){H - z)ip + + j^) ^ . 

The functional integral is over supervectors ip and ■(/', whose components are commuting and 
anticommuting (or Grassmann) fields |6). Thus, in addition to the (physical) particle-hole 
space, an auxiliary "boson-fermion" (bf) space is introduced. To average a product of n Green 
functions, we further enlarge the field space by forming the tensor product with C" {i.e. we take 
n copies of the field space). We do not distinguish between "advanced" and "retarded" fields, 
as we shall make use of the aforementioned symmetry relating with G^ . 



5 



To take advantage of the particle-hole symmetry of the BdG Hamiltonian, we find it conve- 
nient to tensor the field space with yet another factor C^. This extra "spin" degree of freedom, 
referred to by the name of "charge conjugation" (cc) space, is introduced by rearranging the 
quadratic form of the generating functional as follows: 

= i){H-z)'4} + ^'^{-C-^HC-z)i,^ 
The superscript T here denotes the supertransposition operation. If we define 

*=7i(40' »^;^(^ -*'-^-'). 

we have 

i){H - z)il) = ^{H -w)^ , 

where w = S3 2;. We use the notation for the Pauli matrices acting in the charge conju- 
gation space and ctj for the ones acting in the particle-hole space. As is easily checked, the two 
supervectors \1/ obey the symmetry relations 

where 7 is defined by 

j=lph® {Ebb Si - EpF ® iSa) ® 1„ 

with Ebb (Eff) being the projector on the bosonic (fermionic) space. Note that 7 and C 
commute. The trick of doubling to will help to make the symmetries of the nonlinear a 
model manifest at all stages. 

To summarize, the generating functional for averages of products of Green functions can be 
written as 

Z[J] = (^j V^V^ exp j (fx {i^{H - w)* + * J J*) ^ (1) 

with sources J, J that obey the same symmetry relations as \E', 4". The structure of the operators 
in this expression is the following: 

H = Hph (S) 16/ <?> Icc ® Iji , 

w = lp/,0l6/(g)S3 0diag(2;i,...,2;„) 

where zj = Ej — irjj arc complex numbers with negative imaginary parts {—rjj < 0). In using 
this tensor product notation, we will usually omit trivial factors of unity. 

After this digression about notation, we go about deriving the effective field theory for the 
disordered superconductor. First of all, decomposing H into regular and stochastic parts as 
H = 1-10 + V{x)ai, we perform the average over disorder in the generating functional Z: 

Z[Q] = I V^V^ I VVexp I d^x ( z\l/(Ho - w)* + i*CT3*y(a;) ^ 



J V^V^ J VVexp J (fx ^z*(Ho -w)* + i*CT3*y( 
J D*D*exp J (fx ^i*(Wo -w)*- y(*c73*) 



2v^ 

2 
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The next step is to make a Hubbard-Stratonovich transformation, which allows to extract 
the relevant "slow modes" of the problem. We are interested in keeping the diffusive modes, 
i.e. those modes that arise from two-particle channels undergoing multiple scattering with the 
exchange of momenta smaller than the inverse of the elastic mean- free path, i = vpT. Isolating 
these modes is a standard procedure |^ which is conveniently performed in Fourier space and 
has been included in Appendix H for completeness. It is shown there that the contribution from 
the slow modes to —v^{'^a3'^)^2 can be written as 

STrC(-<z)C(<z) 

kl<9o 

where go = 1/^7 and 

k 

One might think that the magnetic term in the Hamiltonian Tio suppresses the slow modes that 
couple to the vector potential A. It will turn out, however, that the magnetic field of a vortex 
is not strong enough to cause complete suppression of these modes. Therefore, we will keep all 
the slow modes appearing in the expression for C{q). 

To enable the treatment of the quartic term, one now introduces a Hubbard-Stratonovich like 
decoupling. This procedure reduces the quartic terms to bilinears, at the expense of introducing 
additional fields Q: 

exp |^-z;2^STrC(g)C(-9)^ = 

J VQexpY, (^STr Q{q)Q{~q) - ^STr Q{q)C{~q) 

The symmetry relation 

STrQ^-^cra = STr cra^^^-^Q^ = STr cr3(7-iC-i^')(-^7 C)g^ 

= -STr C7 Q^7"V3C"^** = STr aij Q'^j-^ai'i'^'cTj , 

is accounted for by subjecting the supcrmatrix Q to the linear condition 

Q = ai7 Q^7"'^^i ■ (2) 

After integrating out the quasiparticle fields ^ and switching back to the coordinate repre- 
sentation, we obtain 

Z[0] = JvQexp-S[Q], 

S[Q] = -J (fx (^^STr Q2 _ isTr ln(Ho - + ia^Q/^T)^ . (3) 

The integration domain for ^' and Q has to be chosen with care, in order to ensure the conver- 
gence of all integrals. A detailed discussion how to choose integration manifolds correctly can 
be found in [ p^ . 

In summary, the problem of computing disorder averaged Green functions has been reduced 
to considering the Euclidean field theory with effective action S[Q]. The sign difference compared 
to Efetov's expression for the effective action comes from the different definition of the supertrace. 
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3 Saddle point approximations 



The next step in deriving the nonhnear a model is to find the saddle points of the effective 
action (|^), and to incorporate the fluctuations around these. Unfortunately, because the gap 
function A varies in space, the full saddle point equation is too complicated to solve directly. 
However, in the dirty limit tA ^ 1, the minimum of the action functional can be found by a 
two-step procedure jll) . In that case the scales set by the disorder and by the superconducting 
order parameter are well separated, so that one can perform two minimizations in sequence. 

To outline this section, the strategy adopted from is as follows. At first, one neglects 
the gap function A and the deviation w of the energy from the Fermi level. One varies the 
resulting effective action and finds the corresponding saddle point manifold. Then, fluctuations 
inside this manifold are considered; they couple to the gap function and to the energies w. The 
resulting low-energy effective action is varied once again inside the first (high-energy) saddle 
point manifold. The partial differential equation obeyed by the optimal configuration turns 
out to coincide with the equation derived by Usadel for the Gorkov Green function of a 
superconductor in the dirty limit. 

3.1 First saddle point 

Let us write the BdG Hamiltonian for a clean field-free superconductor as 



A=0 



/loCTs + A , A = (ReA)cri - (ImA)cr2 , 



where /iq = — V^/2m — In the absence of the gap and at the Fermi energy, variation of the 
action functional S[Q\ gives the following saddle point equation: 

Q{x) = —{x\{ho + iw3^3 + iQ/2T)-^\x) , (4) 

where we kept the positive infinitesimal rj, which is crucial for distinguishing between the physical 
and unphysical solutions. We assume /i ^ l/r. Then, since {x\{hQ + irj)^^\x) is understood as 
—inv, equation is formally solved by a diagonal matrix Q — diag((7i, (72, ■■■), with qi = ±1. 
To choose the signs correctly, we note that the expression on the right-hand side of the saddle 
point equation relates to the Green function for the disordered system in the self-consistent Born 
approximation. The disorder preserves the causal (i.e. retarded versus advanced) character of 
the Green function, and therefore the sign of qi must coincide with the sign of the imaginary 
part of the energy. This singles out the solution 
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3.2 Fluctuations around the first saddle point 

In the limit A ^ 0, A — > and w 0, the action (^) is invariant under the transformation 
Q{x) I— > TQ{x)T^^ , with T any supermatrix that is constant in space and consistent with 
the symmetry condition (|^) for Q. Obviously, the saddle point equation must have the same 
invariance, and its solution will therefore be degenerate. Indeed, the saddle points lie on a 
manifold that can be parametrized by Q = TQqT~^, satisfying = 1. For /it ^ 1 the low- 
energy excitations are given by slow transverse fluctuations preserving the constraint Q^{x) = 
Qq = 1. To incorporate these, we put Q{x) := T{x)QqT^^{x) where T(x) is determined only up 
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to "gauge" transformations T{x) i— > T{x)h{x) that leave Qo fixed: h{x)Qoh~^{x) = Qq. When 
taken modulo such gauge transformations, the matrix T{x) varies slowly with x. Fluctuations 
that are not of this form are massive and will hence be neglected. 

We now insert Q{x) = T{x)QaT~^{x) into the expression for S[Q] and expand with 
respect to A and w, and up to second order in gradients of Q{x), neglecting higher-order deriva- 
tives. In this way the following low-energy effective action governing the transverse fluctuations 
of Q is obtained: 

S[Q] = j (fx STr {D{VAQf - 4»(A + uj)a^Q) . (5) 

Here, D = vj^r/S is the diffusion constant of the normal metal and Va = V — ieA[a3, ■]. The 
details of this calculation can be found in Appendix ^ 

It will turn out to be of importance that, for w = (and any A and A), the action functional 
S[Q] is invariant under transformations 

Q{x) ^ TQ{x)T-\ if T = lph<»t (6) 

and t is constant in space and obeys t — 7(i^^)^7^^. The latter condition means that t runs 
through an orthosymplectic Lie supergroup OSp(2n|2n). 



3.3 Second saddle point: Usadel equation 

We now want the optimal configuration of the superfield Q for a nonvanishing order parameter 
A and magnetic vector potential A. To find it, we require S[Q] to be stationary with respect to 
variations of Q that preserve the constraint = 1 . Such variations can be parametrized by 

SQ{x)^e[Xix),Q{x)] , 

where X satisfies X = —crijX'^j'^ai so as to preserve the symmetry (|^). The stationarity 
condition SS = then yields 

^VA(QVAQ) + i[Q,(A + c^)a3] = 0. (7) 

With proper interpretation given to Q, this coincides with Usadel's equation for the Gorkov 
Green function - a fact first understood in jl^ . By a slight abuse of terminology we shall refer 
to (0) also as the Usadel equation. 

Our plan now is as follows. In the remainder of the current subsection, we are going to 
solve the Usadel equation for a vortex with cylindrical symmetry (at a; = 0). To do that, we 
will introduce cylindrical coordinates (r, tp, z) and look for a z-independent solution Qo {r, if) . 
Because of the invariance (^, this solution is degenerate and, once again, we are dealing with 
a manifold of saddle points TQo{r,ip)T~^, generated by the action of T e OSp(2n|2n). In the 
next section, we will argue that the low-energy excitations of that manifold are given by Q{x) — 
T(z)QQ(r, ip)T~^(z), with T{z) = Iph t{z). The fluctuations in the r and (p directions are 
strongly suppressed by the coupling to the pairing gap of the superconductor and by transverse 
quantization. By inserting the last form of Q{x) into the effective action (^ we will obtain yet 
another effective action, this time for the matrix T{z). 

The first step is to solve the Usadel equation at u; = 0. This was done before in and 
we redo it here in a way tailored to our purposes. Recall, first of all, the invariance of the 
Usadel equation (at w = 0) under orthosymplectic transformations Q{x) i-^ TQ{x)T^^^ with 



9 



T = Iph ^ t and t E OSp(2n|2n). Now, far from the vortex, any solution with finite action S[Q] 
must become covariantly constant (VaQ 0) and adjust to the pairing gap A = | A|crie~"^°'^ 
of the superconductor in the bulk {[Q, Acts] 0). These conditions in conjunction with ~ 1 
determine Q at infinity to be Qoo — cr2e^"^°'^ , up to multiplication by a sign. Note that Qoo is 
stable under symmetry transformations: Qoo = (Iph ® t)Qoo(lp/i ® t)~^ . On the other hand, on 
approaching the vortex core, the unique value of Qoo is expected to expand into a manifold of 
solutions of the Usadel equation, as a result of the orthosymplectic invariance. To remove the 
resulting degeneracy, we hold the solution on the axis of the vortex (r = 0) fixed at the metallic 
saddle point Qo = (73 (8) S3. Given the two limiting forms at r = and r = 00, it is natural to 
look for a solution of the form 

Q(r, if) = g{r)a3 ® E3 + fir)a2e-''^^^ ® l^c , (8) 

which interpolates between Qo and Qoo- Note that this ansatz is compatible with the symmetry 
condition (||). To fulfill the constraint Q^ = 1, the functions g{r) and f{r) are required to satisfy 

f{r)+g'ir) = l. 

Note also that /(r) must vanish linearly at r = in order for Q{r, Lp) to be nonsingular there. 

Next, we insert our ansatz into the Usadel equation, thereby obtaining a differential equation 
for g{r) and /(r). As is shown in Appendix]^, this equation reads 

Ur{rfdrg~rgdrf)+{^^-2eA^^ fg^^\A\g. (9) 

To complete the formulation of the problem, we need to specify the gap function |A(r)| and 
the (/3-component of the magnetic vector potential, Atp{r). In principle these quantities are to 
be calculated self-consistently [p^ from the function /(r). It is common knowledge, however, 
that their essential features are as follows, (i) The gap function vanishes on the symmetry 
axis of the vortex, increases linearly with r for small r, and saturates at the bulk value of the 
superconductor over a scale given by the coherence length ^. (ii) The magnetic field carried by 
the vortex integrates to one fluxon h/2e and is localized in a region of linear size equal to the 
penetration depth A. 

For simplicity, we do not insist on an exact self-consistent solution but adopt the following 
closed-form approximate expressions ||l^ (for an exact treatment, see Rcf. pll ): 

„ , , , 1 Ki(r/X) , , , , r/i 

which capture the features listed above. The parameter a is a numerical constant whose value 
will be specified later on. Notice that the vector potential A is everywhere nonsingular, as the 
Bessel function Ki{r/\) behaves as A/r near r = to cancel the pole 1/r. Note also that, given 
the value of the gap Aq far from the vortex, the coherence length ^ no longer is a parameter at 
our disposal. Rather, ^ has to be determined self-consistently: because A is given by / through 
the BCS gap equation, the length scale ^ fed into the differential equation (^) must agree with 
the characteristic scale of the solution emerging from that equation. 

To proceed, we solve the constraint -I- = 1 by setting / — sin and g — cos 0. If we also 
introduce the dimensionless variable u = r/f , the equation for 9{u) takes the form 

u cos a 



u du \ du J \ K J 2 D ^/u' 
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with K — A/^. Let us now specialize to the case of an extreme type-II superconductor, where 
K ^ 1 and the Bessel function can be approximated by Ki(u/ k)/ k « 1/u. In this Hmit, all 
explicit length scales are removed from the equation - and the self-consistency requirement of 
equal input and output length scales is satisfied - by setting 2^^Aq/D — const. Without loss of 
generality we take the numerical constant to be unity, so that 

e = V^/2Ao . 

(Note that a different choice of numerical constant would only amount to a redefinition of the 
parameter a.) This result says that disorder reduces the "clean" value of the coherence length 
Co = vf/ttAo to a "dirty" value ^ ^ (^^o)^^^- The Usadel equation then simplifies to 

„ 9' sin 20 ucos6' 
~9" + —r- = , . 10 

As was stated before, the boundary conditions are cos6'(oo) ~ .g(oo) — and sin6'(0) = /(O) — 0. 
Without loss of generality we put 6(0) ~ 0. 

For large values of u, we can get an approximate analytical solution as follows. The boundary 
condition at infinity admits the possible values 6{oo) = (2m + l)7r/2 with m an integer. Writing 
0{u) = (2m + l)7r/2 + h{u) and approximating sink by h in the asymptotic regime, we get 

h" + h'/u={-l)"'h. 

The solutions of this second-order differential equation are combinations of Bessel functions. For 
even to, the solution that goes to zero at infinity behaves asymptotically as 

h{u) - u-i/^e-" , 

while for odd to, the solutions have the asymptotic form 

h{u) ~ u^^^^ cos(u + uq) . 

When inserted into the functional (|^), the latter solutions can be shown to produce a radially 
integrated Lagrangian density 



g {r)rdrdif> ^ dr 







which diverges linearly. On the other hand, the solutions for even m, which decay exponentially, 
have a finite radial action. Hence the "good" boundary values at infinity are 6{oo) — tt/2 + 27rn 
with integer n. This conclusion is confirmed by a second argument: in Appendix ^ we show 
that only the solutions with tan ^6{oo) > can be analytically continued to the normal metal 
solution, which has to emerge at energies large compared to the gap Aq. 

Since these solutions with n G Z are local extrema (and, in fact, minima) of the action, one 
might naively think that all of them have to be taken into account, according to the principles of 
saddle point evaluation of integrals. However, the solutions with n 7^ do not appear anywhere 
in the large body of literature on the Usadel equation and its applications. What's the reason? 
The answer is that in the conventional kinetic theory the Usadel equation is based on, one argues 
that the solutions with n 7^ are in conflict with causality. Indeed, g(r) has an interpretation 
[ p7| as the local density of states (LDoS) at radius r, normalized to the density of states of the 
normal metal. By the properties of the cosine, the solution for n = is expected to yield a 
function g = cos 9 that is everywhere positive, but for n 7^ regions of negative value inevitably 
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appear. In other words, the n ^ solutions fail to meet the physical requirement of positivity 
of the density of states, or causality of the Green function. 

If this argument does not seem totally convincing - on the grounds that only the total (not 
the partial) LDoS is required to be positive - there exists also a functional integral reason why 
the n 7^ solutions must be discarded. A careful look reveals that for d{oo) = 7r/2 + 2?™ the 
field Q winds n times around the singularity of the logarithm in equation (|^). Consequently, 
the solutions with n ^ are not smoothly connected to any choice of integration manifold ||l^ 
for Q (before saddle point approximation) that is permitted by the required convergence of all 
integrals. Put differently, the solutions with n ^ arc not accessible by deformation of the 
integration manifold without moving across singularities of the integrand. This uniquely singles 
out the solution with n — 0. In summary, we adopt the two boundary conditions 

61(0) = , and ^(oo) = tt/2 . 

We have yet to fix the value of the parameter a in equation (^0|). To do that, we follow Ref. 

and choose a in such a way as to achieve exact self-consistency far from the vortex. Standard 
theory tells us that the order parameter A = X{ipiip^) (where —A is the BCS coupling constant) 
is obtained by summing the thermal Green function Jp^(-0|(T)?/'f (0))e*"'^(iT over all Matsubara 
frequencies cui = {21 + 1)T/-k. At temperature T = 1/13 = the sum condenses to an integral 
over the (imaginary) energy axis, and by expressing the thermal Green function through the 
function /(w, E) — sni6{u^ E) one obtains 

/•oo poo 

A(m)/Ao = / sin 6 {u, ~iuj)duj / / smd{oo, ~iuj)duj . (11) 

The function sin 9{u,E) is determined by the Usadel equation at energy E (see Appendix |C|): 

^„ 6' sin 26* wcose* E 

-0" \ — = , - z— sm 6* . 

u 2u^ Vm^ + ^0 

The prescription of Ref. now is to expand the solution around the asymptotic limit 



(oo, —iuj) — Aq/^ Al+uj^ , cos 6{oo, —ioj) — uj / \J A^ + oj"^ , 



find the correction of order 1/u^, insert the corrected solution into equation ([ll|), and equate 
the result for A(w) with the ansatz A(u)/Ao = uj \Ju^ + d?- = 1 — a^/2u^ + ... made above. This 
leads to a linear condition for a giving = 7 r/2 . 



We have solved the differential equation (10) with a = n/2 numerically, using Maple. Fol- 



lowing again Ref. |9|], we integrated the equation from u = to u = oo, with the derivative 9'{0) 
so adjusted as to hit the desired value 7r/2 at infinity. The solution is shown in Figure 1. In the 
next section we will see how this solution enters into the expression for the coupling constant of 
the nonlinear a model. 



4 The one-dimensional nonlinear a model 

We are in the process of deriving a low-energy effective theory describing diffusive quasiparticle 
motion in a disordered vortex. The final step, which we now embark on, begins by taking the 
solution of the Usadel equation and perturbing it by fluctuations that vary slowly in space. 
Recalling T(r, (p, z) = \ph ® t{r, ip, z), we put 

Q = TQvseuieiT-^ = gaa ^ t^^t'^ + /aae"*'^"^ 1 . 
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e(u) 




Figure 1: The solution 9{u) of the differential equation ( p^ ) with the boundary conditions 
9{0) = and 9{oo) = tt/2. The value of the parameter was chosen to be tt/2. 



This expression is now inserted into the effective action (^. Let us first evaluate the "kinetic 
term" STr (WaQ)^- After some calculation we find that the magnetic field disappears completely, 
the corresponding terms involving either a trace over traceless matrices in the particle-hole space 
or a supertrace over a supermatrix proportional to the identity. We thus obtain 

STr (VaQ)^ = STr (VQ)^ . 

The explicit form of the gradient is 

VQ - {Vg)a3t^3t-^ + gas^itEst-^) + V {f cj2e-'^^') . 

By the cyclic invariance of the supertrace and STr 1 = 0, the kinetic term reduces to 

STr(VQ)' = 5'STrp,,xfc/xccx«(lph<»V(iI]3t"i))' 
= 252STrb/xccxn(V(iI]3t-i))' . 

In the last expression the particle-hole space has disappeared, its only remnant being the factor 
of 2 multiplying the restricted supertrace. By introducing a reduced Q-field, 

q = ti:^t-^ , 

we can write more simply 

STr [VQf = 2.g2 STr {Vqf . 
The other term in S[Q] works out to be 

STr (A + uj)a3Q = STr Ao-aQusadd + 9 STr tj(lph §5 iSat^-^) = 2g STr ujq . 

The effective action for the field q then reads 

S[q] ^-^Jd'x STr {g^Di^qf - Aigcuq) , 
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To identify the low-energy modes, we use cylindrical coordinates and write the operator 
controlling the fluctuations around a constant background g-field as 

-.9 V.gV = -7^ irj—r 27T T~ r-9 (^)^ 

oz'^ OLp^ rg'^yr) or \ or 

This operator is self-adjoint in the i^-space with measure g^{r)rdrdipdz. The spectrum of the 
first term on the right-hand side is continuous, while the second and third term have a discrete 
spectrum. [The latter is seen by making a similarity transformation 

—— (rg^{r) — ] ^ -r'^/^ g-\r)—rg^ {r)—r-^/^ g-\r) = -— + V{r) , 

rg^{r) dr \ dr J dr dr ' dr"^ ' 

and showing that the effective potential V{r) is confining.] Thus the transverse, i.e. r and if, 
excitations of the g-field are gapped (with the gap being of the order of the transverse Thouless 
energy 13/^^), and the relevant low-energy physics is captured by taking g to be a function of z 
only. Doing so we arrive at the following low-energy effective action for q = q{z): 

S[q] = / dz STr {C2D{d,qf - 4^ClWg) , (12) 

where the dimensionless coefficients Ci and C2 are given by radial integrals, 

C„ = 2 / cos" d{u) udu (n = 1, 2) , 

which are of order unity. To calculate their precise values one needs the function cos ^(u) which 
is determined by self-consistency. If we take 6{u) to be the solution of equation (|lO|) with 

= tt/2, we get Ci = 3.16 and C2 — 1.20. This completes the derivation of the low-energy 
effective theory. What we have obtained is a one-dimensional field theory that belongs to the 
general class of nonlinear a models. The expression (|l^) for its action functional is the first 
important result of the present paper. 

To put the result into context, we elaborate on the definition of the field space. Recall 
that the field q{z) is a 4n x An supermatrix q = iSat"^ with t{z) e OSp(2n|2n). Clearly, the 
expression for q is invariant under right multiplication of t(z) by a matrix h{z) that fixes E3 
(/iSa/i"^ — S3). By inspection, one easily sees that the group of such matrices is isomorphic 
to GL(n|n). Thus, the field space of our nonlinear a model is a quotient OSp(2rt|2n)/GL(n|n). 
Let us emphasize that the present theory is distinct from the nonlinear a models that were 
introduced by Efetov Q in his work on disordered metallic wires. In fact, what we are dealing 
with here is one of the new "normal-superconducting" universality classes (namely class C) 
identified in Ref. To demonstrate the distinction clearly, let us take a close-up look at the 
field space. 

Since 7 = {Ebb ® Si — Epp §5 iT,2) ® In gives E3 = —7(23)^7"-^, the Lie group condition 
t = 7(i~^)^7~^ translates into q = —jq'^'j^^. A useful way of representing q = tE3t~^ is by 

(1 + ZZ)(1 - -2Z{1- ZZ)-^ 

2Z{1-ZZ)-^ -{1 + ZZ){1 - ZZ)-^ 

where Z and Z are complex 2n x 2n supermatrices. The condition q = —7q'^J~^ is equivalent 
to 



ZFB ; - [-{Z^^f {Z^P f 
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We see in particular that Z^^ , Z^^ are symmetric matrices, while Z^-^, Z^^ are skew. More- 
over, stability of the functional integral with action (^2|) requires 

Z^^^-(Z^^)t, i^^=+(Z^^)t, l-(Z^S)tz^«>0. 

When complemented by the OSp(27i|2n) action induced on Z, Z by g ^ tqt^^, the above 
properties mean that Z^^ and Z^^ parametrize certain symmetric domains that are well-known 
in mathematical physics. Technically speaking, the _Bi?-sector is isomorphic to the noncompact 
symmetric space S0*(2n)/U(ri) (called type Dili in Cartan's notation), while the Fi^-sector is 
isomorphic to the compact symmetric space Sp(2n)/U(n) (called type CI). In the terminology 
of Ip^ , these facts are succinctly summarized by saying that the field space of our nonlinear a 
model is a Riemannian symmetric superspace of type ZJIIIICI. In contrast, Efetov's models are 
defined over the Riemannian symmetric superspaces of type AIII|AIII, BDl\Cll, and Cll\BDl. 

Before turning to the analysis of the one-dimensional field theory, let us briefly discuss its 
zero-dimensional limit. Such a limit is obtained if the vortex is taken to be isolated, with 
Neumann boundary conditions {dzq){0) — {dzq){L) = at the two end points z — Q and z — L. 
(The boundary conditions on the microscopic quasiparticle wave functions at an insulating 
boundary are of course Dirichlet, but the reduction to the effective field theory turns them into 
Neumann.) At small energies (jJ, the dominant contribution to the functional integral comes 
from the "zero mode" q{z) = qo (independent of z). Spatially varying field configurations have 
a Lagrangian density that scales with the length of the vortex as L~^. By comparing terms in 
the action functional, one sees that they are negligible if 

Lu < D/L^ . 

This condition, which delineates the zero-dimensional (or ergodic) regime, is familiar from dif- 
fusive metals: ergodic behavior sets in at times much longer than the diffusion time L^/D or, 
equivalently, at energies much smaller than the Thouless energy Etyi — D/L^. 

In the zero-dimensional limit, the functional integral over q{z) reduces to a definite integral 
over the zero mode qo: 

Z[J] = j Dqo exp (—iTrSTr Li;go/2(5 -|- sources) , 
1/5 = . 

To obtain any specific observable, one keeps track of the coupling to the source terms introduced 
in the starting expression (|l|). For the special case of the 71-level correlation function, the 
corresponding definite integral was computed for all n in . The result is 

, , 1 sin27ri?/(5 , , 

for the mean density of states as a function of the energy difference E from the Fermi level, and 

VI N /sin7r(ei - £2) sin7r(£i £2)^ 

>2 ei,£2) = 7 ^ r^— {e^ = E^ 5) 

V 7r(£i-£2) 7r(£i+£2) / 

for the two- level cluster function I2 ■ These results agree exactly with the expressions predicted 
in P] on the basis of a random matrix hypothesis for systems of symmetry class C . 
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How does this compare with the results expected on the basis of the traditional approach 
due to Eilenberger and Usadel? The Usadel equation by itself, i.e. without the nonperturbative 
extension worked out here, simply yields 1/5 for the mean density of states at energies compa- 
rable to the level spacing. While this answer is correct asymptotically at 3> 5, it ignores the 
oscillations seen in (|l3|) and, in particular, it misses the fact that the density of states vanishes 
at £■ = 0. (Note that the Usadel equation predicts secular variations of the density of states 
over a characteristic scale given by the transverse Thouless energy D /(^^ . For our purposes, this 
is a "high-energy" scale which we shall ignore.) 

Why does random matrix theory predict the fine structure of the density of states correctly, 
whereas the quasiclassical Usadel equation captures only the asymptotics? Perhaps the answer 
is not entirely obvious and the following remark helpful. What is well understood is that random 
matrix theory provides a valid description of disordered systems in the ergodic regime, where 
the physics is governed by the zero modes or constant modes of some field theory. [Recall that 
the ergodic regime is the region of validity of (|l3|).] One might object, however, that a vortex is 
a spatially localized object. The solution of the Usadel equation for a disordered vortex is far 
from being constant but has a characteristic spatial dependence. One might therefore wonder 
why a zero mode or random matrix argument has any chance of being correct in the present 
context. 

Our analysis has already given the answer to this question: although the solution of the 
Usadel equation does vary in space, the symmetry group operating on it is independent of spacel 
As we have seen, the random matrix answer is recovered by integrating over the orbit swept 
out by the (constant) action of the symmetry group. A key feature is the mechanism by which 
the theory avoids the superficial conflict between constancy of the zero mode and locality of 
the vortex solution: the orbit generated by the group action contracts to a single point Qoo = 
TQoaT~^ (thereby making the symmetry ineffective) as we move away from the vortex. 

To conclude this section, let us mention that similar statements apply also to an SNS junction 
with phase difference 0i — 02 = tt (symmetry class CI). In that case, too, the random matrix 
predictions of Q are expected to hold as they stand. 

5 Transmission along the vortex 

Charge transport through a diffusive metallic wire at zero temperature undergoes a crossover 
from Ohmic behavior to strong localization as the length of the wire is increased. A similar 
crossover is expected for the transport properties of the normal quasiparticles of a disordered 
vortex. In the case of metallic wires, the precise form of the crossover functions for the electrical 
conductance was computed from the solution of Efetov's nonlinear a models in one dimension 
in |2^. (Note that an error made in these computations for /3 = 4 was corrected in ]2l|.) 
We will now embark on a similar project for the vortex, by analysing the one-dimensional field 
theory found above. Note that for the purpose of testing for localization of normal quasiparticles 
in a vortex, charge transport is not an experimentally viable option, as the system is short cut 
by dissipationless charge flow through the superconductor. However, we can imagine probing 
the transport properties via a measurement of the heat flow along the vortex. Whether such a 
measurement is experimentally feasible will be discussed in detail in a separate publication. 

In this section we will show that the normal quasiparticles in a disordered vortex are subject 
to a phenomenon called "weak localization" , which is well known from disordered metals with 
time-reversal symmetry and conserved electron spin. Although in normal metals weak local- 
ization is destroyed by the application of a magnetic field, in a disordered vortex there exist 
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additional modes of quantum interference that arise from Andreev scattering of the quasiparti- 
cles and are not suppressed by the presence of a magnetic flux. The zero-dimensional precursor 
of weak localization in systems of symmetry class C was first noticed in a transfer-matrix cal- 
culation by Brouwer and Beenakker ||2^ and was qualitatively explained in p^ . 

As a quantitative measure of quasiparticle transport and localization in a disordered vortex, 
we study the following two-particle Green function: 



Tkk' 



{x, x'; E) = Y^ (4a G+„, {x, x'- E)v]^.2,G'„,^[x', x- E) 



where v is the left-right symmetric velocity operator 

= 2^ (a^ - I - 

The superscripts ± distinguish between retarded and advanced Green functions, and a, a' are 
the particle/hole indices. The angular brackets indicate an average over disorder. 

If -0 = {i^a}a=p,h is a solution of the BdG equations (with some energy E)^ the bilinear 
jk{x) = '^^{'>paVka'4'a){x) expresscs the fc*'' component of the probability current carried by a 
quasiparticle with wavefunction ip. This current is conserved (^i^dkjk = 0), as an immediate 
consequence of the BdG equations. On the same grounds, the tensor Tkk' satisfies the current 
conservation laws 

E ^kk' {x,x'-E) = = Y. ^ ^kk' [x, x'-E). 

k k' 

This suggests that Tkk'{x,x';E) can be interpreted as a correlation function of the conserved 
probability flux carried by quasiparticles with energy E. 

The nonlocal tensor Tkk'{x,x'; E) carries more information than we are actually interested 
in here. To eliminate some of the details, let Pz and Pz' be two planes that intersect the axis of 
the vortex at right angles, and consider the integrated longitudinal correlation function 

t{E) := J (fx J d^x' T33{x,x';E) . 

By current conservation, this is independent of z and z'. In fact, t{E) is a measure of the 
"transmittance" of the vortex for quasiparticles of energy E, and plays a similar role as the 
electrical conductance in a metallic wire. 

To set up the functional integral for t{E), recall the relation 

Ga'ai^\x'-,E) = -^CapCa'P'G'ljij,{x,x';-E) , 

013' 

which follows from the particle- hole symmetry of the BdG Hamiltonian. Introducing two replicas 
and choosing to represent the Green functions in the bosonic sector, we write 

rkk'{x,x';E) = ^ C„,C„,,, ((^^Vf 4^2^^ 

aa'PP' 
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where the second hne introduces the charge conjugation space, by making use of the relations 
established in Section ^. The multiple superscript comprises the boson-fermion, replica and 
charge conjugation indices, in this order. The double angles mean an average over disorder and 
over the configurations of the fields 'i'(x). The symmetry breaking term in the action functional 
is given by 

^ = Iph (8) Ifc/ ® S3 ® Ea^ . 

Now, after Hubbard-Stratonovich transformation, integration over the quasiparticle fields 
and saddle point approximation, we obtain 

where 

a 

For a consistency check on this calculation, note that the Usadel equation yields 

k k 

in agreement with the current conservation law satisfied by the tensor Tkk'{x, x'\ E). 
The last step is to use the parametrization of Section |4| for the matrix Q: 

Q = g{r)a3 ® t{z)^3t-\z) + /(r)a2e-''^"^ , 

and to integrate over the planes Pz and P^' . We then obtain 

riE)=co{.f'''^''{z)f'''^'\z')) 

where j = qdzQ and 



Co 



For a diffusive vortex of short length, L, we expect the transmittance t{E) to obey Ohm's 
law (t 1 / L) . Our goal now is to verify this law and compute the weak localization correction, 
i.e. the next term in an expansion in the ratio of system size over the localization length. To 
that end, we write the action functional (1121) in rescaled form: 



S[q] = I ds STr {{dsqf - 2teXq) 



16L 

s being the dimcnsionless coordinate along the vortex. The ratio L/^ioc, where 

will later be identified as a localization length, plays the role of the coupling constant of the 
one-dimensional field theory. The parameter e, 

2CiEL^ 

e 



C2D 
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represents the quasiparticle energy measured in units of the Thouless energy E'th = D/L'^, and 
A is the matrix 

A = Ifc/ (g) (73 (g) S3 . 

In keeping with what was said in the introduction, the functional integral is supplemented by 
open boundary conditions at the two ends of the vortex: 

g(0) = S3 = q{l) . 

It is well understood that those arc the correct boundary conditions to use when the system is 
in good contact with infinite particle reservoirs. (Roughly speaking, the deviation of q{z) from 
S3 measures how often a typical Feynman path of the quasiparticle visits the location z. A 
diffusive path that arrives at one of the end points s = or .s = 1, most likely exits the vortex 
and disappears into the large phase space of the reservoir, with vanishing probability of return. 
This draining effect by the coupling to the reservoirs causes the fluctuations of q away from S3 
to be suppressed at an open boundary.) On moving the two planes and Pz' to the ends of 
the vortex (z = 0, = L or s = 0, s' = 1), the expression for the transmittance r becomes 

r{E) = - (^^y (j^ii.«22(o)^.B22,sn(i)^ _ 



5.1 Pcirametrization of q 

Recall that the supermatrix q satisfies the orthosymplectic condition 

q = --yg'^'y-i , with 7 = {Ebb Si - Epp ^ I2 • 

The kinetic term STr (dq)"^ in the action functional S[q] is invariant under global transformations 

q{z) 1-^ tq{z)t^^ with t € OSp(4|4). This symmetry is broken by the potential term STr \q at 
finite energy e ^ 0. We will see that increasing the symmetry breaking coupling e causes 
crossover to Efetov's nonlinear a model with unitary symmetry (or (3 = 2), thereby eliminating 
the weak localization correction in the limit e — > cx3. 

The effect of symmetry breaking is best captured by parametrizing q as 



q = toti'Esitot 



where the matrices to, ti decompose into blocks in charge conjugation space as follows: 
The 4x4 blocks Zq, Zo, Z\, Z\ arc chosen to have the following substructure in replica space: 

- i)- --Co : 

Zo = { ~i) , z.-f^"^ 



Here all quantities Z/^, Zy, , z^ and their partners with tilde are 2x2 supermatriccs with 
symmetries that are determined by g = — 7g^7~^. The rationale behind this choice is that to 
commutes with A {tQ^Xto = A), so that the symmetry breaking term reduces to 

STr Xq = STr A<otiS3(toti)"^ = STr Xt^E^t^^ , 
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which simplifies the calculation appreciably. A similar parametrization was used by Altland, 
lida and Efetov |Q in their study of the spectral statistics of small metallic particles subject to 
a symmetry breaking magnetic field. The boundary conditions at s = 0, 1 translate into 

Zo(0) = Zo(0) = Zo(l) = Zo{l) = 0, 

with identical relations in force for Z\^Z\. Using these, we can linearize the current j = qdz<i in 
Zq^ Z\ at the boundaries, and the expression for the transmittance becomes 

r = -(6oc/i)' (5.4'''' (0) 5.4''''(1)) • 

5.2 Expansion of the kinetic and potential terms 

We will now use standard field theoretic perturbation theory to compute the leading terms of 
the Taylor expansion in L/^ioc for t. The first step is to expand the action functional S'[g] up 
to fourth order in the fields. Setting q — toiiS3(toii)^^ we obtain for the kinetic part of the 
action 

STr {d,qf - STr {^^d,t^+tl\t^^d,t^)t^ , S3]' . 

By inserting the explicit expressions for Iq and ti in terms of Zq, Zq etc. and using the fact that 
Zq and Zi are orthogonal to each other with respect to the supertrace, we get for STr (dgq)'^ up 
to second order 

-8STr (dZadZo + dZidZi'j , 
and a fourth order contribution given by 

- SSTr (^dZodZoiZiZi + ZqZo) + dZodZo{ZiZi + ZqZo)) 

+ SSTr (^ZodZo{Zj_dZi-dZiZi) + ZodZo{ZidZi~dZiZi)j (14) 

+ 8STt (^{dZaZi)^ + {dZoZif -dZidZiZiZi-dZidZiZiZi'j . 

For simplicity we have set d = dg. For the potential term we have the expansion 
STr Xq = 4STr (za^a - ^v^v) + 4 (STr (za^a)^ - STr (zvZv)^) + ... , 

which depends only on the matrix entries of Z\^Z\. 

5.3 Propagator and diagrams 

The quadratic part of the field theory action defines the bare propagator. We write 
S\q\ = ^ ds STr {{dqf - 2ie\q) = ^2 + ^4 + • ■ ■ 

where ^ 

<5'2 = tIy^ / As f STr (dZdZ) + ieSTr (za^a — ZyZy)^ , 
2L Jq \ I 

and Z in this expression means Zq + Zi. From the quadratic part of the action it is clear that 
the modes of type Za are massless, while the modes of type and za have an imaginary mass 
term given by e. In the limit e 00, the latter modes will therefore be completely suppressed. 
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Figure 2: Graphs for the propagator {Z'^^Z^ " ) 

We have yet to discuss the symmetries of the matrices Z and Z. It is easy to see by 
hnearization that the symmetry q = —jq^j^^ is equivalent to the condition 



X = -7X^7^1 on X = 



Z 
Z 



On decomposing 7 into blocks in charge conjugation space, 

we see that the 4x4 supermatrices Z = Zq + Zi and Z = Zq + Zi must obey 

Z = -aZ'^ , Z = -Z^a , 

where a = (o'a)^/ (g) I2 is the superparity matrix, with value +1 on bosons and —1 on fermions. 
Recalling the definition of the supertranspose T and denoting the superparity of a by | a |, we 
equivalently have 

^ a| |6| 2"^'^ Z'^^ ( ^1 ^^1^1 ^1^1 1^1 ^^'^^ 

In this notation, the kinetic term reads STr dZdZ = J2abi-^y''^9Z''''dZ'"'. 

The basic Gaussian integral governing the perturbation expansion is given by 

Z-\s)Z''-'{s'))^ JvZVZ Z-''is)Z''-'is') e^v j\sY^i-lpZ,,d^Z,}j 

This formula applies to the massless modes Zq. For the modes Zi {i.e. and Zv) one needs to 
replace —d^ by —9^ + ie and —d^ — ie respectively. The index structure on the right-hand side 
is determined by the super skewsymmetry of the matrices Z and Z. It is seen that there are 
two different ways of propagating the matrix indices, which we represent graphically as shown 
in Figure ^. 

The space dependence of the propagator is given by 



for the massless modes, and by 



s(l - s') for s<s', 
s'(l - s) for s > s' , 



_ 1 J sinh\/Ies sinh\/ze(l — s') for s < s' , 

Viesinh V«e I sinhi/Ies' sinhA/Ie(l - s) for s > s' , 



21 



for those with mass squared equal to ie. Note that these Green functions satisfy Dirichlet 
boundary conditions at the ends s = and s = 1, as required. The Green function Dv(s, s') — 
{—d^ — ie)^^{s, s') is obtained from D/\(s, s') by the substitution e — e. 

We are now ready to calculate the transmittance of the vortex up to one-loop order. What 
we want is the quantity 

r = -teoc/i)' (az„^i^^2(o)azf -^1(1) X (1 - 54))^ , 

where 6*4 is the part of the action functional that is quartic in Z and Z, and the vacuum 
expectation value is taken with respect to the Gaussian action 82- Si contains terms from the 
kinetic and the potential part in S. 

Let us first compute the leading contribution (tree level). Using the expression for the 
massless propagator we obtain 

TO = -teoc/i)^(9Zo^^'^^(0)azf^^^(l))^ 

= -{^ioc/L)dsd,,-Do{s,s') =6oc/i. 

s=0,s' = l 

This is the equivalent of Ohm's law for the electrical conductance. Note that dsds'Do{s, s') is 
in fact independent of s and s', as is expected from conservation of the quasiparticle probability 
flux. The Ohmic limit tq = £,ioc/L gives support to the interpretation of ^loc as a localization 
length. In a diffusive wire one would have ^loc ~ Ni with N the number of channels open at 
the Fermi energy or, equivalently, ^loc ~ I'AD with A the cross sectional area of the wire. In 
the present case we saw that 

This means that the role of the cross sectional area is here taken by ~ ^ ^o^- 

We are now going to compute the weak localization correction, which is of order (i/^ioc)°. 
First, we notice that the "external legs" dZ^ ^' ^(0) and dZ^ ^' ^(1) can only couple to terms 
containing at least one Zq and one Zq. Therefore the fourth order terms coming from the 
potential, as well as the terms on the third line in (^), can be omitted from ^4 for present 
purposes. In Appendix ^ we show that the contributions originating from the terms on the 
second line of (|l^) cancel each other, so we need only consider the terms coming from the first 
line: ^ 

Ti = III (^dZ^^'^'^{0)dZ^'^-^\l) ds STr {ZZdZ^dZ^ + ZZdZ^dZo) 

where we again used the notation Z = Zi) + Zi and the fact that Zq and Zi are orthogonal with 
respect to the supertrace. 

Next, we draw the one-loop diagrams that potentially contribute to the transmittance. The 
ones associated with the term STr ZZdZ^dZ^ are displayed in Figure ^. Those associated with 
STr Z ZdZi^tdZo are similar, except that the loop is placed in the lower part of the graph. The 
first two diagrams, which contain a freely running index, vanish by the usual mechanism of 
supersymmetry: the fermion loops come with a minus sign relative to the boson loops. 

Let us now compute the contribution from Diagram 3. In the course of this computation it 
will become clear that the diagrams with massless loops (type Zq) are zero. First we consider 
the following contraction of the fields: 

^(-l)l"l ldZ^^'^^{0)dZ^^{X)dZ^ dZl'^^Z^'' 

abed \ 
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Figure 3: One- loop diagrams potentially contributing to the transmittance. 



which yields the expression 



E 



dZ„ [0)8 Zq {s)j^(^dZ„ 



{s)dz^'^^\i))Jz^'''' 



{s)Z'''^\s) 



The propagator associated with the loop, {Z^^''^{s)Z'^'-^^{s))q, consists of two terms, one con- 
taining a factor [—l^'^lS-^^'-^^S'^''' (Diagram 1) and another one with —S^^''^5'^'^^ (Diagram 3). 
The first term vanishes by supersymmetric cancellation, = 0, as stated above. (Note 

in particular that all contributions from Zo-loops cancel by this mechanism.) The second gives 



si=0 



dsBoisi, s) Da(s, s) dsds 



Do(s,S2) 



Diagram 4 gives a very similar contribution, the only difference being that D\/(s,s) replaces 
D/\(s,s). It should also be clear now that if the external legs are contracted with the fields Z 
and Z, the loop is of Zq type and yields a vanishing contribution. 

The same result is obtained when the diagrams arising from STr dZodZoZ Z are considered. 
Although there might in principle be extra contributions from the expansion of the functional 
integration measure, we show in Appendix ^ that this is not the case. 

In summary, we obtain the following result for the weak localization correction to the trans- 
mittance: 







ds 5si5sDo(si,s) (Da(s,s) +Dv(s,s)) dsds^Bois, S2) 



si=0 



S2 = l 



^ ( sinh\/ies sinh\/ie(l ~ _|_ sinh \/—ies sinh\/— ie(l — s) 
Jo I \/ie sinh\/ie y/—ie sinh ^/~ie 



1 / coth Vie coth \J —ie 
'2 1^^^+ ^ 
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Using the expansion of the hyperbohc cotangent around x = 0, 

1 X 
coth X — — h — + . . . , 
X 3 

we find that the weak localization correction for quasiparticles of zero energy is ri(e = 0) = —1/3. 
In the limit e ^ oo the correction vanishes as l/v^- The crossover between the two regimes 
takes place at e ~ 1, corresponding to energies of the order of the Thouless energy Erh ^ D/L'^. 

To finish this section, let us briefly recall the physical picture explaining the origin of the 
weak localization correction computed above. Weak localization in normal-/superconducting 
systems with a magnetic field arises from modes of interference that are present already in 
the single-particle Green function. The relevant diffusive processes are described by single- 
particle Feynman paths containing a closed loop which is traversed twice. (This is the physical 
interpretation of the third and fourth diagram in Figure |^.) To make such a process phase 
coherent in the presence of a magnetic field, the charge state of the quasiparticle during the 
second traversal must be opposite to its state during the first traversal (the charge is reversed 
by Andreev reflection). In that case, the total magnetic phase accumulated in the loop vanishes, 
and the process survives disorder averaging. The net effect of such processes is to suppress the 
density of states and reduce the transmittance of the vortex (weak localization). 

The preceding discussion assumes the energy of the quasiparticle to be zero or negligibly 
small. At finite energy E, the phase coherence between the two charge states in the closed loop 
is reduced by an additional dynamical phase, which is proportional to E multiplied by the time 
typically spent in the loop. For diffusive dynamics, this time is estimated to be of the order of 
L^/D in an open system of size L. Therefore, phase decoherence due to i? > becomes effective 
at energies of the order of the Thouless energy D/L^. 

In the limit E 3> D/L'^, states of opposite charge can no longer interfere, so the problem re- 
duces to considering BdG particles only, or BdG holes only. The surviving modes of interference 
are the same as those of a metallic wire subject to a magnetic field. It is well known that the 
latter system has a vanishing weak localization correction, and one must go to two-loop order 
to see the onset of localization. The field theory describing such systems is Efetov's nonlinear a 
model with unitary symmetry {(3 = 2). What one expects, then, is that the field space of q{z) 
becomes isomorphic to Efetov's model space (/? = 2) when the modes Zi, which are killed by 
the limit E —^ oo, are set to zero. Such an isomorphism in fact exists, and will be constructed 
in the next section. 

6 Strong localization 

To introduce the subject of this section, we briefly recall how the conductance of a disordered 
quasi-one-dimensional normal metal varies with the length L of the system [Q. Provided that 
L is not too small, the disorder averaged conductance g{L) depends only on a single scaling 
variable: 

g{L) - F^(L/6oc) , 

where ^loc is a characteristic scale for the onset of localization. The scaling function is 
specific to each Wigner-Dyson universality class (with symmetry index /3 = 1, 2, or 4). Let the 
normalization conventions for ^loc be such that Fp has the Ohmic (or small x) limit Fp{x) = 
1/a; + ... in all cases. Then, to next order in the small-x expansion, 

Fp{x) = --cp + 0{x) , 
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where the constant cp, the weak locahzation correction, is cp = 2/3, 0, —1/3 for /3 = 1, 2, 
4. Looking back at the result of the previous section, we see that weak locahzation for the 
low-energy quasiparticles of a disordered vortex is qualitatively the same as for metallic wires 
with symmetry index /? = 1. 

In the opposite limit of large x the scaling functions fall off exponentially: 

F^(x)"-^A^exp(-x//3), 

which is an expression of the phenomenon of strong localization that occurs in quasi-onc- 
dimensional systems. Similar behavior is expected for the transmittance r of the vortex. The 
purpose of the present section is to confirm this expectation and work out the precise value of 
the exponent governing the exponential decay law. 

Recall the formula expressing r as a functional integral over fields Z: 



For the purpose of extracting the largc-L asymptotics, it is best to abandon the functional 
integral formulation and pass to an equivalent quantum mechanical description. To do that 
we interpret the dimensionless coordinate s as an imaginary time. Since the Lagrangian of our 
system is 

(6oc/16L)STr {-{dsqf + 2ieXq) , 
the evolution with "time" [i.e. length) will be governed by the Hamiltonian 

W = --^A+^STrAg, 

where A is the Laplace-Beltrami operator of the field space G2/H2 = OSp(4|4)/GL(2|2) with 
metric tensor — STr {dq)^/8. The open boundary conditions on the fields mean that we want a 
quantum transition amplitude from Z = 0ats = 0toZ = 0ats = l. Moreover, by the rule 
that X for a Euclidean path integral with Lagrangian density + V{x) translates into the 

difi'erential operator m~^d/dx in the quantum theory, the expression for r becomes 



T=(Z = 



j^Bl,B2 Hr)B2,Bl 

Uq e Uq 



z = o 



where D^^-^^ = d/dZ^^'^^ and i)^^"^^ = d/dZ^^^^'^. 

Supersymmetry guarantees the ground state of H to be exactly at zero. By inserting a 
complete set of cigcnstatcs of H, we sec that the exponential decay of t for large L is determined 
by the smallest eigenvalue of Ti. such that the corresponding eigenfunction ip has nonzero matrix 
element {0\Dq^'^^\iP). Unfortunately, computing the spectrum (or even the lowest nonzero 
eigenvalue) of Ti for all e is a rather difficult task. The only cases we can do easily are e = and 
e = 00. These are the cases we will restrict ourselves to here. 

Consider the case e = first. The Hamiltonian we are facing is minus the Laplacian on 
G2I Because the Laplacian commutes with the action of the symmetry group G2, its spec- 
trum is in principle computable by the powerful machinery of Lie group theory. There exists, 
however, a trick we can use to further simplify the problem. Although we are forced to introduce 
two replicas to average a product of one retarded and one advanced Green function when the 
quasiparticle energy E is nonzero, in the limit = = e we can make do with just a single 
replica. The idea is to use again the BdG particle- hole symmetry to express G~{E) by G~^{—E), 
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and then to represent one Green function in the bosonic sector and the other one in the fermionic 
sector. This can be done without breaking the supersymmetry of the formahsm if and only if 
= 0. Developing the rest of the theory in the same manner as before, we arrive at a reduced 
description in terms of a "small" q field taking values in G\jli\ = OSp(2|2)/GL(l|l). (The 
external derivatives Dq^'^^ and Dq^'^^ in the transition amplitude formula for r are replaced 
by odd derivatives D^^ and D^^ .) 

To obtain the lowest nonzero eigenvalue of minus the Laplacian on Gi/Hi with a minimal 
amount of mathematical effort, we adopt the pedestrian approach of choosing a suitable coordi- 
nate system and expressing the Laplacian as a differential operator in those specific coordinates. 
A particularly useful coordinate system is provided by the following polar decomposition: 

q = hai^sa h , a = exp ( ^ 1 , ^ = I q -0 j , 

Here 0, tp are real commuting variables, and C, are complex anticommuting. We refer to 6 as 
the "radial" coordinate. The metric tensor in this coordinate system takes the form 

-STr {dqf/8 = d9^ + sin^ 26 {d(p + CdC/2i + CdC/2i)^ + sin^ 6 dCdC . 

Now, there exists a standard formula for the Laplacian on a Riemannian (super)manifold: 

ij 

if the metric tensor is '^C'ffy (O'^C"' , and g{^) = SDet(gij(^)). Given this formula one can 
write down a lengthy expression for the Laplacian on Gi/Hi. However, there is yet another 
simplification and we do not need to go through the trouble of writing the full expression. 
The point is that we can construct the desired eigenfunction of A by repeatedly acting with 
symmetry generators on a corresponding radial eigenfunction f{9). It is therefore sufficient to 
study the spectrum of the radial Laplacian, which from the above expression for the metric is 



A# 



sin26'\ d /sin26'\ d 



sin^e dO \sin^e d9 



The low lying spectrum of this operator can be constructed very simply by trial and error. 
Making the substitution x = cos 26 we have 

_A*/4 = -(l-.)A(l+.,-i, 

and it is easy to verify that this has the eigenfunctions 

/o = l, fi = l-x, f2 = {l-x){l + 3x) , 

with eigenvalues Aq = 0, Ai = 1, and A2 = 4. The pattern we see suggests that the complete 
spectrum of —A"^ /A might be with n an integer. (This guess is in fact confirmed by more 
informed Lie algebraic considerations.) In any case, it is clear that the lowest nonzero eigen- 
value is Ai = 1. By applying to /i = 1 — a; an odd generator of the symmetry group Gi, 
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one can explicitly construct an odd eigenfunction tp of — A/4 with eigenvalue +1 and matrix 
element (0|-D^^|V') 7^ 0. Since the Hamiltonian is H = — LA/2^ioc, we predict the decay of the 
transmittance to be 

t{E = 0) - exp (-2L/6oc) 

for quasiparticles with zero energy. 

The other limit E 00 can be done more quickly, by drawing on known results. As we saw, 
this limit eliminates the modes Zi, and what remains are the modes Z^, which are recalled to 
have the matrix representation 

By the off-diagonal structure of these matrices, the field space decomposes into two separate 
pieces, which are isomorphic copies of each other by the orthosymplectic relations Zq = —ctZq 
and Zf) = —ZqCt. We can write the Lagrangian in terms of one of these copies, say Zq 
and its partner Zq , and account for the other one by multiplying by 2. The kinetic term 
-(6oc/16i) STr then takes the form 

(6oc/i) STr (1 - z,z+)-'d~z^{l - z+S^)-'dz+ . 

Recall that Zq and Zq are complex 2x2 supermatrices unconstrained by any symmetry relations 
that do not involve complex conjugation. This means that the object at hand is the Lagrangian 
of Efetov's nonlinear a model for metallic wires of the Wigner-Dyson symmetry class /3 = 2, in 
the specific paramctrization given by 

and with a renormalized value of the localization length. This model has been much studied in 
the literature ^5, 19, Taking the result from there, we find 



t{E = 00) - exp (-L/4aoc) . 

Thus the localization length at e = is predicted to be shorter than at e 3> 1 by a large factor of 
8. This strong energy dependence already announced itself in the weak localization correction, 
and now we see that it becomes even more dramatic in the region of strong localization. 



7 Concluding remarks 

In this paper we have established a field theoretic description of the diffusive motion of normal 
quasiparticles along the core of a vortex in a dirty, extreme type-II, s-wave superconductor. 
The technique used was an extension, due to Ref. |pl| , of the supersymmetric method devel- 
oped by Efetov and others for the purpose of studying localization and mesoscopic phenomena 
in normal metals. A comprehensive account of this extension and its application to normal- 
/ superconducting mesoscopic systems will appear shortly |2^ . 

In a simplified picture, one may view a vortex core as a normal conducting wire surrounded 
by superconducting walls. The quasiparticles in such a system undergo Andreev scattering 
and partake in interference processes that are not present in a disordered normal metal. This 
gives rise to different spectral correlations and transport properties. For a short, isolated vortex 
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the correlations of the quasiparticle spectrum (in the universal low-energy limit) have been 
verified to agree with those of a random matrix of symmetry class C . For an open system and 
quasiparticle energies below i?Th = D jL? , the "transmittance" defined in Section |^ exhibits a 
weak localization correction. Strong localization sets in for vortices of length L > ^loc- The 
localization length ^loc for a quasiparticle with energy above the Thouless energy E^h is eight 
times larger than the localization length of a low-energy quasiparticle. 

These results were obtained by reducing the effective field theory for the quasiparticles to 
a one-dimensional nonlinear a model. In the process, we had to solve a saddle point equation, 
which turned out to be the Usadel equation. The saddle point approximation is valid if ^ 1/t 
and 1/r 3> Aq, the latter being known as the dirty limit. The reduction to one dimension 
(the degree of freedom along the vortex) can be done for quasiparticle energies smaller than the 
transverse Thouless energy -D/^^. 

It is a well known fact that nonmagnetic, isotropic disorder has no significant effect on s-wave 
superconductivity (Anderson's theorem). This ceases to be true when the order parameter is not 
spatially homogeneous. Then, a self-consistent field treatment is required. The self-consistency 
problem for the disordered vortex has been studied both in the framework of Landau-Ginzburg 



theory and by using the Usadel equation |10|. Both approaches show that disorder with 



an elastic mean-free path i changes the coherence length from its clean value to a dirty value 
C - (^^o)^/^- We have seen that this scaling of the dirty coherence length comes about quite 
naturally, by using no more than dimensional analysis and some minimal assumptions about 
the shape of the condensate. 

We did not implement an exact self-consistent treatment here. Instead, we took the order 
parameter profile A(r)/Ao from a simple one-parameter family of functions and chose the free 
parameter so as to achieve self-consistency in the region of approach to the asymptotics far 
from the vortex. We believe that an exact treatment will lead only to minor corrections in the 
numerical constants Ci and C2 that enter the coupling constants of the nonlinear a model. In 
any event, the qualitative picture that emerges from our results, namely, a crossover between the 
low-energy regime with weak localization present in the "transmittance" along the vortex, and 
the high-energy regime without weak localization, is a robust feature that will not be modified 
by a more quantitative treatment. 

For typical values of a superconductor with a normal metallic diffusion coefficient of the 
order of 20cm^/sec, and a critical temperature around lOK, a rough estimate of the localization 
length gives a range 1 — 10mm. This corresponds to an elastic mean-free path of the order of 
5nm and a coherence length of about lOOnm. For a sample of 1 — lO/im thickness, the Thouless 
energy of the vortex will be in the mK range. Because of the large localization length, the weak 
localization correction is a relatively small effect. Still, the predicted weak localization correction 
to the transmittance along the vortex should be visible as a low-temperature anomaly in the 
heat conductivity of an ensemble of vortices between two electron reservoirs. We shall discuss 
this effect in detail in a forthcoming communication. 

Note added. After completion of this manuscript, we learned of related work by Skvortsov, 
Kravtsov, and Fcigel'man [ pTj . 

Acknowledgment. We are much indebted to Alexander Altland for valuable discussions 
and input at many stages of the present work, and for locating Ref. ||^. The initial spark for this 
research came from Gerd Bergmann, who a number of years ago asked one of us (M.R.Z.) how 
the disordered vortex relates to the conventional classification scheme for disordered electron 
systems. 
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A Decoupling of the slow modes 

Consider the quartic term 

which appears on averaging over the random potential V{x). We want to isolate in it the 
collective modes with small momentum, \q\ < qo, which are to be decoupled by a Hubbard- 
Stratonovich transformation. To do that, we Fourier transform to the momentum representation: 

Jd^x {^{x)as^{x))' = J ^^^^(fc0a3*(W(fc3)a3*(-fci-fc2-fc3). 

There exist three independent ways of pairing two fast single-particle momenta to form a slow 
two-particle momentum q: 

*(A:i) *(A:2) *(A:3) - ^2 - A:3) 

a) k —k + q k' —k' — q 

h) k -k' -q -k + q k' 

c) k k' —k' — q —k + q 

Term o) can be decoupled trivially, producing no more than energy shifts that can be absorbed 
by a redefinition of the chemical potential. The other two terms can be rearranged in the 
following way. For term b) we have 

^ *(/c)a3*(-A;'-g)*(-/c + g)a3*(fc') 

k^k' yq 

= J2 ^{k)(T3^{-k' - q)^^{k')a3^^{-k + q) 

k,k' ,q 

= ^{k)a3^{-k' - q) (-*(fc')7"'C) (73 (7C-i*(-A; + q)) 

k,k' ,q 

= ESTr|^^*(-fc'-(?)*(A;')a3j ||X^M/(-fc + g)*(fc)a3j . 
Here we have used Cu3C~^ = —03 and the relations 

which follow from ^E* = C^^^^^ and 'I' = — 'I'-^C^^7 (sec the main text). The term c) is the 
"nice" one and is easily put in the very same form, by using only the cyclic invariance of the 
supertrace. Thus, the quartic interaction J S'x (^0-3^)^ is replaced by 

J2 STr a-QKil) , 
\Q\<qo 

where ( is given by a sum of dyadic products of the fields ^ and ^: 

C(g) = ^*(-fc + g)*(fc)a3. 
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B Effective action 



We want to extract from the functional S defined in (||) a "low-energy" effective action for slowly 
varying fields Q{x) — T{x)QoT^^{x) generated from the saddle point Qq ^ (J3 (E) S3. To do 
that, we insert Q{x) — T(x)QqT^^{x) into S: 



S[TQoT-^] = ^ I cPx STr In ( T^o - ^ + '-^TQoT-^ ) {x,x) , 



2t 

and expand with respect to uj, A, and A, and up to two gradients of Q{x), which are assumed 
to be "small" quantities. 

Let us first extract the gradients. For that, we set w, A and A to zero and consider 

5o = i y d'x STr In (^a3(-VV2m - ^l) ^ '^TQoT-"^^ {x, x) . 

This cannot readily be expanded as it stands, as T(x) is ambiguous by multiplication on the 
right by an arbitrary field h{x) with h{x)QQh~^{x) = Qq. To fix the ambiguity, we write T[x) 
as 

T{x) = (7e^(=") where X{x)Qo + QqX{x) = . 

The matrix g is chosen in such a way that gQog~^ coincides with the local "mean" of Q{x), while 
X(x) parametrizes small deviations from that mean. The space dependence of g is assumed to 
be so weak that it can be neglected for the present purposes. 

When Q = ge^ QQe~^ is inserted into So, the constant matrix g drops out by the cyclic 
invariance of the supertrace, and expansion of the logarithm in powers of the small quantity 
X{x) gives 



So 



i / (fx STr In 



"2^ - + ^e^Qoe"-^ ] {x, x) 



1 

'8^ 



(fx I (fx' STr 



2m 



M+^Oo) {x,x'){X{x')-X{x)) 



M+^Qo) {x',x){X{x)-X{x')) + 



We assume the elastic mean- free path £ to be much larger than the Fermi wave length [kp^ ^ 1). 
The Green functions then are well approximated by 



2m 



m exp(^z/ci?|a; — x'\ 



-\x-x'\l1l 



27r 



|a; — x'l 



and Taylor expansion of the fields 
leads to the integral 



fc=i 



d^x' 



„'|2 



2m 



2t 



{x,x') 
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Assembling factors we obtain 



TO- 

On = 



Unr- 



■ J Sx STr {\/Xf 



By recalling the formulas v = mkp/2'K'^ — m^^/27r^T for the local density of states and D = 
P' I'iT for the diffusion constant in three dimensions, we can write this as 

So^^J d^x STr (VX)2 . 
To put the result in invariant form, we use 

STr (VQ)2 ^ STr (V(ge^Qoe-^g-i))' = STr (Ve2^)(Ve-2^) = -4STr (VX)^ + ... . 
This gives 

So = -^Jd^x STr (VQ)2 . 

which is the kinetic part of the effective action (|^) . 

The dependence of the Lagrangian on the magnetic field is easily restored by the following 
symmetry argument. Consider a gauge transformation 

e 

The transformation law induced on the composite field Q ^ 'i'^a^ is 
so the gradient VQ transforms as 

VQ ^VQ + i{W(p)[a3,Q] . 
This determines the covariant derivative Va acting on Q to be 

Va = V - ieA[a3, ■] . 

Therefore, the unique way of making the effective action gauge invariant is to substitute 

STr (VQ)2 ^ STr {VaQ)^ , 

which restores the dependence on the magnetic field without any calculation. 
Finally, we expand to extract the terms linear in A and cu. These are 

Si = l f d^xSTr (z^-^,+ ±q] {x,x)a3{A-uj) . 
2 J \ 2m 2t J 

By using the saddle point equation (— V^/2m — /x + «(5/2r)^^(a;, x) = —mvQix) and the relation 
CT3A — — A(T3 we immediately get 



S'l = ^ y Sx STr Q(A + w)ct3 



(Note that, since Q and A behave in the same way under gauge transformations, the term 
STr QAcTa is gauge invariant, as is STr Qujuy,.) Summing the terms S'o and Si we arrive at the 
effective action (||). The sign differences compared to come from: a conjugation of Q by 0-3 
(which changes the sign of the off-diagonal terms of Q in the particle-hole space) , and an overall 
change of sign of Q. 
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C Differential equation 

To solve the Usadel equation for a vortex configuration, we made the ansatz Qq = g(r)a3 (8) 
E3 + /(r)tT2e~*'^'^^ ® Icc in Section |37^ . The purpose of this appendix is to derive the differential 
equation obeyed by the coefficient functions / and g. 
To begin, we recall the Usadel equation for Q: 

DVa(QVaQ) = ~i[Q, (A + uj)a^] . 

Since all of the quantities Qq, Aa^, A = |A|(Ti exp(— iyjua) and uj = ET,^ are diagonal in charge 
conjugation space, we may specialize to the sector where E3 equals unity and work with the 
simpler ansatz 

Q = 9ir)a3 + /(r)(T2e-'^'^^ . 

(Alternatively, we may imagine absorbing E3 into g.) The right-hand side of the Usadel equation 
then becomes 



i[ga3, AfTa] -^[/a2e-^'^"^^a3] - {i\A\g + Ef) 2aie- 



The evaluation of the left-hand side of the Usadel equation is more laborious. From = 
V — ieA[a3, ■] we have 

Va(QVaQ) - V(QVQ) - ieV{Q[Aa3, Q]) - ieA[a3, QVQ] - e^A'icTs^Qias, Q]] . 

We now express this in cylindrical coordinates {r,ip,z). For the term with two derivatives we 
get 

V(gVQ) - ^drirQdrQ) + ^d^iQd^Q) + d,iQd,Q) 

= -drirfdrf + rgdrg) -f -drirfd^g - rgdrf)cj^e-''^-' 
r r 

+ ((.9^3 + /a2e-^'^'^^)/e^^-^ai) 



-dr{rfdrg - rgdrf) + ^ ] iaie 



The term fdrf + gdrg vanishes hy + g^ — 1. Taking the magnetic vector potential to be 
purely azimuthal, A = eipAip{r), we obtain for the part linear in A 

IP IP 
- leV {Q[Aa3,Q]) - ie[Aa3,QVQ] = --d^{Q[A^cj3,Q]) ~ -[A^a3,Qd^Q] 

Op A TP A 

= -^^^ ii9<^3 + /a2e-^^-^)/e''^'^Vi) - ^ h, {gas + /a2e-^'^'^^)/e'^-Vi] 
r 

The term quadratic in the vector potential gives 

-e^A^ [(T3, Q[a3, Q]] = -e^Alia^QasQ - Qa^Qas) = 46^^2/5 iaie"*'^"^ . 
Gathering all terms and dividing by icrie"*''"^^ we arrive at the equation 

^drirfdrg - rgdrf) + (^J - 2e^^^ fg=^ (|A|g - iEf) , 
whose solution for = is discussed in the text. 
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Figure 4: Definition of tlie angle cj) for a) cj^ — 77^ — > and b) — 77^ — < 0. 

D Choosing between solutions 

The full saddle point equation (including A and oj) is difficult to solve exactly if the gap function 
varies in space. However, for the case of a uniform gap a solution is easy to come by. This solution 
is constant and describes the region far from the vortex. For values of the quasiparticle energy 
OJ large compared to Aq, the solution should approach the normal metal solution. This fact 
will now be used to analytically continue the metallic solution (w/Aq S> 1) to the zero-energy 
solution for a superconductor (w = 0). 

If we ignore the charge conjugation degree of freedom, which is not essential here, the full 
saddle point equation reads: 



In the absence of any spatial dependence in the gap A, this is solved by Fourier transforming 
to momentum space and taking Q to be constant: 



Here, w > and 77 > 0. We choose the gauge where A = Aui and diagonalize the matrix 
M = 7Acr2 — (a; — 777)(T3 by a similarity transformation with a matrix T: 



The sign of the first eigenvalue, A, of M can be chosen so that A —oj + irj when w/A 00. 
We then have: 



where the angle (p varies from vr to 7r/2 as goes from infinity to (see Figure The imaginary 
part of A remains positive for any (positive) value of to. 



The matrix that diagonalizes M is T = (1 + &cti)/\/T~~62^ .^^j^^ 5 = + \J (^^) - 1- 



It is important that the imaginary part of b is negative for any value of w/A. Note in particular 
that for w = we have h = —i. The saddle point equation now takes the form: 
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The physical solution is given by Q = TQT ^ = sign(ImA) x 0-3 = (T3. Therefore 

Q = ^ ^-3^ = + . 

Unphysical solutions (which break causality) also exist. 



By comparing to the ansatz used in Section 3.3 to solve the Usadel equation, Q — a^i cos9 



(72 sin 6*, we conclude that b = —i tan |. In particular, for = where b = —i, we get tan | = 1- 
Therefore, 9 = ■|(mod27r) is the zero-energy solution for the homogeneous superconductor, in 



agreement with the asymptotics found in Section 3.3 from a different argument 



E Vanishing vertices 

In this appendix we show that the terms on the second line in equation ^ do not contribute to 
the transmittance t{E) in one-loop order. Let us consider the first and second term, which give 



^aZo^''^^(0)aZo^^'^'(l) STr ZodZoiZidZi - dZiZi) 

As before, only the contractions of type 1 and 2 in Fig. ^ are nonzero, and we obtain 

Do(6, s) (D/,(s, s) - D;,(s, s)) Do(s, i) + (A ^ V) , 

where the dot over an argument represents a derivative acting on it. Because the first derivatives 
of the propagators Do(s, s'), D/\(s, s') and Dv(s, s') are discontinuous at s = s', these terms are 
ambiguous and have to be defined by some point-splitting procedure, 6.17. 

DA(i,s) = lim-(D;,(s + -,s)-DA(s--,s) = ^ ^ ' 



a^o a V 2 2 / 2sinh v«e 

Since the propagators are symmetric in their arguments, we have 

Da,v(s,s) = D/\^v(s,s) , 

which establishes our claim. In the same way it is seen that the one-loop contributions from the 
third and fourth term in equation ^ cancel each other. 

F Invariant integration measure 

The field space of the nonlinear a model we have derived is a coset space G/H with G = 
OSp(2n|2n) and H = GL(n|n). The space G/H is realized by matrices 

where g € G, and S3 is fixed by the elements h of the subgroup H: 

It is easily seen that the symmetry relation q = —jq^'j^^ corresponds to the condition 
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which defines the group G. Let Dii{q) denote the superintegration measure (on G/H) invariant 
under q ^ g^qg^^ for all go £ G and normalized by / D^{ci) = 1. Our goal in this appendix is 



to express D^{q) in the parametrization of q introduced in Section 5.1 



Quite generally, if we are given a parametrization with (super)coordinates and metric 
tensor 

STr(dg)2=^d^g.,(C)de^ 

then the invariant integration measure is locally (i.e. modulo boundary ambiguities) expressed 
by 

D^Ji{q)=Di SDeti/2(5,,(0) , 

where is a (suitably normalized) Euclidean measure. Doing the computation this way requires 
knowledge of all matrix elements of the metric tensor and is more complicated than is necessary. 
Since we are parametrizing q in the form 

g(0 - .9(0^33-^0 , 

a more efficient approach is to calculate the measure directly from the Jacobian of the map 

{i\---,e''')^i^gm^g-\0 ■ 

However, this is still not the best method as it involves the group element g twice, once on the 
left and once on the right (by the inverse). 

The most efficient approach is to turn to the "coset space picture" . Given some basis 
Xi, . . . ,Xdim of = Tff{G/H), we assemble the coordinate functions ^* into a matrix-valued 
function X = J^i ■ G/H V and consider some map 

X^g{X)H, 

which sends X into the coset space G/H . Its differential is 



[XAX)^- 
as 



giX + s- dX)H . 



Because the integration measure on G/H we wish to express by X, is invariant under left 
translations, we may multiply by g{X)~^ and switch to a map Tx : V ^ V defined by 



TxidX) = - 
ds 



g{X)-'g{X + s-dX)H 



s=0 



g{Xy^giX + s-dX) 



s=0 



The subscript V means projection of a Lie algebra element A e Lic(G') = Lic(i?) + V onto 
A-p G V. Let J{X) = SDetTx be the Berezinian (or super jacobian) of the transformation. 
Then the invariant measure DjjL^q) is simply expressed by 

DXJ{X) . 

We now wish to apply this general formalism to the parametrization 



1 



(a = 0,1) 
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This is not possible immediately, as the matrices do not obey the law g = 'j{g~^)'^'j~^ defining 
the group G = OSp(2n|2n) (although q = tQtiT,^ti^tQ^ does satisfy q = — 79"^7~^). To be on 
safe mathematical ground, it is advisable to start from a parametrization 

such that the Ta do satisfy the group law. This is achieved by choosing 



I - z^z^)-^/^ ZM - z^z^)-'/^ 

1 Zc.\ f h+ 

Za 1 M h- 



where /i+ = (1 - Z„Z„)-V2 and /i" = (1 - Z„Z„)-V2. 
The general formalism now tells us 

Dn{q) = D{Zo, Zo) D{Zi,Zi) J{Zo, Zo, Zi, Zi) , 

where D{Zq, Zq) and D{Zi,Zi) are flat Berezin measures, J = SDet T is the superjacobian, and 

T{dZo,...) = {r,-\^'d{T,T,))^ = (T,-'dT,)^+(T,-\T^'dT,)T,)^ . 

Recall the decomposition Z ^ Zq + Zi which, in terms of spaces, we write P = 'Po + 'Pi- Because 
of {Ti^dTi )-p = (rf^dTi)^^ and the general relation 

SDet ( 5 =SDet^^ ^ 



z y V z 



we have SDet T = SDet T' with 

T'idZo, ...) = {n'dT,)^^ + (rr^(To-Mro)Ti)^^ 

If X is any element of Lie(G) we introduce the notation 



X 



(A{X) B{X)\ 
[ciX) D{X)J 



where ^ Lie(il), and e "P. A straightforward calculation then gives 

S(TfMri) = h^dZ^h^ 
C{T^'^dTi) = h^dZ^ht , and 

B ((rf Vq-Mto Ti )p J = h+ (h+dZo /lo - Z, h^dZo h+Z, ) /if , 

c((^r^'ro"^droTi)p„) = h^(hQdZoh+ - Z^h+dZf^h^Z^^h+ . 

By writing out the matrices in explicit form, one easily sees that the Berezinian of the transfor- 
mation dZi h^dZ^h^, dZ^ ^ h^dZ^ht is given by SDet(/i^) SDet(/i]") = SDet(l - Z^ZJ, 
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while the Berezinian of the transformation dZg i-^ hgdZgh^, dZ^ i— )■ K^dZ^h^ is simply unity 
as a result of cancellations due to supersymmetry. 

For similar reasons, the Berezinian of the transformation 

dZp {dZQ — Z^dZf^Z^ jh^ , 

dZg 1-^ h^{dZf^ — Z^dZ^Zi)}!^ 

turns out to be unity. Thus, by collecting factors the invariant integration measure D^i{q) in 
the parametrization q = TqTiT,^Ti^Tq^ is given by 

D{Zo, Zo) D{Zi, Zi) SBetil - Z^Zi) . (15) 

Finally, we switch to the desired parametrization 

It is seen that this amounts to transforming Ti t~> h^T^ h^^ or, equivalently, 

Z^ IiqZj^Iiq, Z^ Hq Z-j^h^ . 

Once again, the superjacobian of this transformation is unity, so that Dj.i{q) in the desired 
parametrization has the form (^5|) as it stands. 

In one-loop order of the perturbation expansion for the one-dimensional field theory, the 
factor SDet(l — Z\Z\) contributes a two- vertex J ds STr ZiZi where a is the regularization 
scale. This vertex does not couple to the external legs dZ^^'^"^ and dZ^^'^^. Therefore we may 
pretend the functional integration measure to be flat, which is what we did. 
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